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SIMULATING  THE  OPENING  OF  FUSIBLE-LINK-ACTUATED  FIRE  VENTS 


Leonard  Y.  Cooper 


ABSTRACT 

Model  equations,  suitable  for  general  use  in  compartment  fire  models,  are  developed  to  simulate  the 
thermal  response  to  arbitrary  fire  environments  of  fusible-link  activated  automatic  fire  vents.  The 
method  of  analysis,  which  focused  on  a prototype  three-element  link-mount  design,  can  be  extended 
to  arbitrary  multiple-element  link  mounts.  Also,  the  equations  for  the  prototype  three-element  link 
mount  are  shown  to  include,  and  be  directly  applicable  to  the  problem  of  simulating  the  response  of 
single-element  simply-supported  links. 

A method  was  developed  to  determine  the  values  of  the  set  of  parameters  that  characterize  a 
particular  multiple  link-mount  design.  This  involves:  1)  the  measured  time-dependent  thermal 
response  of  the  link-mount  design  to  exposure  in  a plunge  test,  the  type  of  test  used  to  characterize 
the  thermal  response  of  sprinkler  links;  and  2)  an  analytic  means  of  determining  values  of  the  design 
parameters  that  yield  an  optimum  fit  between  a solution  to  the  model  equations  and  the  temperature 
data.  The  method  takes  account  of  the  pre-fusing  solder-melting  phase  of  fusible-link  response.  The 
latter  was  shown  to  be  very  important  for  the  link  design  used  in  this  study. 

The  method  was  carried  out  successfully  for  four  plunge  tests  involving  three  link-mount  systems 
and  one  or  two  system/flow  orientations.  In  each  case,  the  model  equations,  with  the  determined 
values  of  the  system  parameters,  lead  to  excellent  simulations  of  the  plunge  test  temperature  data. 
The  resulting  system  parameters  and  model  equations  were  used  successfully  to  simulate  previously- 
reported  time  of  link  fusing  and  vent  activation  in  a real  fire  environment. 


Keywords:  building  fires,  compartment  fires;  computer  models;  fire  models;  mathematical 

models;  vents. 


Automatic  Vent  Opening  Devices:  Fusible  links  and  Drop-Out  Vents 

Ceiling-deployed  fire  vents  are  one  traditional  means  of  venting  smoke  during  fires.  Although  many 
possible  automatic  or  manual  strategies  for  the  opening  of  vents  are  possible,  a common  strategy 
involves  automatic  actuation  by  use  of  a thermally  sensitive  device,  deployed  near  the  ceiling  of  the 
facility,  in  the  immediate  vicinity  of  the  vent  itself.  Two  typical  devices  are  offered  by  vent 
manufacturers,  the  fusible  link  and  the  “device”  of  the  drop-out  vent. 

In  the  case  of  the  fusible  link,  a spring-loaded  mechanism  is  designed  to  open  the  vent  when  restraint 
hardware,  which  includes  the  fusible  link,  is  ruptured.  In  the  event  of  a fire,  the  high  temperature 
smoke/fire  gases  rise  to  the  ceiling  of  the  protected  facility,  turn  and  flow  across  the  ceiling  and  past 
the  near-ceiling-deployed  link,  and  heat  it  by  forced  convection.  If  the  link  temperature  rises  to  its 
fuse  temperature,  TFUSE,  the  solder  that  joins  two  bonded  parts  of  the  link  melts.  As  a result,  the  link 
and,  therefore,  the  vent-opening  restraint  hardware  ruptures,  and  the  vent  opens. 

The  drop-out-type  vent  “device”  involves  a vent  frame,  mounted  in  the  ceiling  of  the  protected 
facility,  whose  opening  is  covered  by  a vent  “dome.”  The  dome  is  formed  from  a plastic  sheet,  and 
is  designed  to  soften,  shrink,  and  drop  out  of  the  frame  when  heated  to  a sufficiently  high 
temperature.  The  dropping-away  action  leads  to  the  open-vent  configuration.  As  in  the  case  of  the 
fusible  link,  the  heating  of  the  vent  dome  is  primarily  by  convective  heat  transfer  from  the  high 
temperature  smoke/fire  gases  that  flow  across  the  ceiling  and  across  the  lower  surface  of  the  vent- 
dome  plastic.  This  work  does  not  address  further  the  subject  of  drop-out  vent  actuation. 

The  Objective  of  this  Work 

The  focus  of  this  work  is  on  vent  system  designs  that  depend  on  fusible  links  for  their  actuation.  The 
objective  is  to  develop  a general  mathematical  model  to  simulate  the  action  of  fusible-link-actuated 
ceiling  vents,  where  the  model  is  suitable  for  general  use  in  compartment  fire  models. 

Simulating  the  Thermal  Response  of  Simply  Supported  Fusible  Links 

Figure  1 illustrates  a typical  fusible  link  design  used  in  fire  vents.  This  involves  two  pieces  of  steel, 
each  with  characteristic  dimension  of  a few  cm,  joined  by  a solder  connection  over  an  area  of  a few 
cm2.  The  link  is  rated  by  the  manufacturer  according  to  the  temperature  at  which  the  solder  melts, 
Tfuse-  Commonly  used  TFUSE  ratings  are  74  °C  and  141  °C. 

As  seen  in  Figure  1,  there  is  a hole  at  each  end  of  the  link.  The  ends  of  steel  cabling  are  typically 
passed  through  these,  and  the  cabling  and  link  become  an  integral  part  of  the  above-mentioned  vent- 
opening restraint.  This  kind  of  link-mount  design  will  be  referred  to  here  as  the  simply-supported 
link  design. 

The  traditional  means  of  predicting  the  thermal  response  of  a simply-supported  link  is  relatively 
straight  forward.  This  is  the  case  since  convective  heating  of  the  link  has  been  typically  assumed 
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to  dominate  the  heat  transfer  process.  In  particular,  the  insignificant  thermal  contact  with  the  cable 
connections  at  each  end  of  the  link  leads  to  negligible  conductive  link-to-cable  heat  transfer.  Also, 
it  is  generally  assumed  that  radiation  heat  transfer  does  not  play  a significant  role. 

The  traditional  method  by  which  the  thermal  response  of  simply-supported  links  are  simulated  is  the 
same  as  the  original  method  used  for  the  thermal  response  of  fusible-link-actuated  sprinklers 
[ 1 ] [2] [3],  where  conduction  losses  from  the  link  to  the  base  of  the  sprinkler  were  neglected.  As  in 
the  case  of  sprinklers,  the  measured  value  of  the  link’s  Response  Time  Index  (RTI)  [4]  and  TpuSE  are 
required  to  implement  the  method.  With  known  values  of  time-dependent  gas  flow  speed  and 
temperature  local  to  the  link,  u(t)  and  Tc(t),  respectively,  the  temperature  of  the  link,  TLINK(t), 
initially  at  the  ambient  temperature,  T^g,  can  be  obtained  from  the  solution  to 

dTLINK/dt  = (Tg  - Tlink)u  /RTI;  TLINK(t  = 0)  = T^g  ( 1 ) 

up  to  the  time,  tp^,  when  T = TrosE.  In  fire  modeling  applications,  Eq.  (1)  is  often  used  to  simulate 
the  response  of  both  sprinkler  and  vent  links  (see,  e.g.,  [5] [6]). 

Vent  Designs  Involving  Complex  Link  Mounts 

Automatic  vent  designs  can  incorporate  the  action  of  Figure  1-type  fusible  links  that  are  not  simply- 
supported,  and  that  can  not  be  expected  to  be  simulated  by  Eq.  (1).  This  is  so,  for  example,  in  the 
case  of  the  vent  design  used  in  large-scale  tests  conducted  by  Underwriters  Laboratories  (UL)  for 
the  National  Fire  Protection  Research  Foundation  (NFPRF)  [7]  [8]  [9],  to  investigate  the  interaction 
of  sprinklers,  vents,  and  draft  curtains  in  industrial  fires,  and  for  Industrial  Risk  Insurers  and  Conspec 
Systems,  Inc.  [10],  to  investigate  the  relative  sensitivity  of  roof  mounted  heat  and  smoke  vents.  As 
will  be  presented  below,  the  latter  vent  is  designed  with  Figure  1-type  fusible  link  with  one  end 
simply  supported  and  the  other  end  welded  and  bolted  to  the  steel  frame  of  the  vent. 

An  important  aspect  of  the  NFPRF-sponsored  tests  was  the  application  by  NIST  of  a computer  fire 
model  that  uses  a computational  fluid  dynamics  methodology  based  on  Large  Eddy  Simulation  (LES) 
analysis  [9].  A goal  of  the  present  work  is  the  development  of  mathematical  models  of  fire  vent 
actuation  that  can  be  readily  adapted  by  a LES  model. 

A typical  fire  scenario  in  the  NIST/NFPRF  Research  Program,  including  a sketch  of  the  vents  used 
(1.2  m x 2.4  m),  is  depicted  in  Figure  2.  As  can  be  seen  in  this  figure  and  in  Figure  3,  the  link  is  part 
of  a three-element  assembly,  around  which  flows  the  elevated  temperature  gases  of  the  fire-plume- 
driven  ceiling  jet.  The  assembly  includes  a 1.2  m-long  structural  steel  “L”-bar  (0.076  m x 0.102  m 
x 0.0048  m thick)  that  spans  the  vent  frame,  a steel  tab  welded  to  the  center  of  the  bar,  and  the  link 
itself,  bolted  to  the  tab.  Above  the  link  and  connected  to  it  is  a cable  connected  to  the  vent  opening 
mechanism.  In  the  usual  way,  fusion  of  the  link  leads  to  opening  of  the  vent. 

In  this  vent  design,  heat  transfer  to  each  of  the  elements,  the  bar,  the  tab,  and  the  link  itself,  of  the 
three-element  link-mount  assembly  plays  a significant  role  in  the  thermal  response  of  the  vent- 
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actuating  fusible  link.  The  next  section  develops  a set  of  model  equations  to  simulate  the  thermal 
response  of  the  three-element  mount.  This  would  be  used  to  predict  fusing  of  the  link  and  actuation 
of  the  vent  for  arbitrary  fire  scenarios.  The  equation  set  is  designed  to  be  versatile,  in  the  sense  that 
it  can  also  be  used  to  simulate  the  thermal  response  of  fusible  links  implemented  in  other  relatively- 
complicated  link-mount  configurations,  and  even  in  the  simply-supported  configuration. 

Simulating  the  Thermal  Response  of  the  Three-Element  Link-Mount  Assembly 

General  comments.  Corresponding  to  the  three-elements  link-mount  assembly  of  Figure  2,  the 
mass,  effective  specific  heat,  and  temperature  of  the  bar,  the  tab,  and  the  link,  are  designated  as  m3, 
c3,  T3.  m2,  c2,  T2;  and  m1?  c1?  T1?  respectively.  The  elements  are  assumed  to  be  physically  small 
enough  and  close  enough  to  each  other  that,  in  estimating  the  convective  heat  transfer,  the 
instantaneous  free-field  gas  flow  local  to  each  of  them  can  be  assumed  to  be  spatially  uniform,  at  the 
same  temperature  and  speed,  TG  and  u,  respectively. 

Convective  heat  transfer  to  the  elements.  Relative  to  the  gas  flow  over  element  i,  the  Reynolds 
number  is 


Re,  ~uLconv,/vair(Tg-i)  (2) 

where  u and  TG  are  the  speed  and  temperature  of  the  gas,  L^oNv.i  is  a characteristic  length 
of  element  i relative  to  estimates  of  gas-to-element  convective  heat  transfer,  and  vAIR(TG.I)  is  the 
kinematic  viscosity  of  air  at  temperature 

Tg.!  * (Tg  + T,)/2  (3) 

The  Rej’s  for  fire  scenarios  and  element  geometries  of  interest  here  are  expected  to  fall  within  the 
approximate  range  100  to  5,000.  Using  a result  for  the  average  rate  of  heat  transfer  to  a circular 
cylinder  in  a crossflow  of  air,  it  is  reasonable  estimate  the  individual  average  transfer  coefficients, 
h,  from  (p.  242  of)  [11]  where 

Nuj  = hjLcQNv  j/k  air(T  g.j)  = 0.48(Rej)/2  (4) 

kA]R(TG_i)  is  the  thermal  conductivity  of  air.  Then,  using  Eqs.  (2)  and  (4)  leads  to 

hi  = ^•'^^(U^CONV,i)1/2kAIR(TG-iyVAIR  ^C^G-l)  (5) 

As  can  be  seen  in  Table  1 

for  20°C  < Tgi<  500  °C: 

(6) 

kAK(TG.I)/vMR  ,/!(Tc.,)  = 6.16(1  ± 0.03)  Ws1,2/m2K 
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Table  1. 

Some  properties  of  air  at  standard  pressure  [12]. 

Tg., 

[°C] 

kA,R(TG.,)  v.\ir*Tg.()  kAIR(TG.I)/vArR  (TG-I) 

[W/(mK)]  [m2/s]  [Ws“/(m2K)] 

20 

0.0251 

15.7xl06 

6.33 

100 

0.0307 

23.6xl0'6 

6.32 

200 

0.0370 

35.5xl0'6 

6.21 

300 

0.0429 

49.2xl0'6 

6.12 

400 

0.0485 

64.6xl0'6 

6.03 

500 

0.0540 

Sl.OxlO'6 

6.00 

Therefore,  from  Eq.  (5),  hj  can  be  reasonably  approximated  as 

hs  = 0.48u'/!LCONVi'1,2kAIR(TG.i)/vA]R'/XTG.I)  = 2.96(u/LCONVj),/’  [Wsl,2/(m2K)]  (7) 

Radiative  heat  transfer  to  the  elements.  In  existing  fire  models,  the  effect  of  radiation  heating  on 
the  thermal  response  of  fusible  links  in  fire  environments  is  not  included.  This  is  true  in  the  case  of 
links  used  to  control  both  vent  or  sprinkler  actuation.  Yet,  in  a proper  analysis  of  link  thermal 
response  the  effect  of  radiation  may  not  be  negligible.  For  example,  consider  the  three-element  link- 
mount  of  Figure  2 where  the  width  of  the  link  is  LC0NV3  = 0.019  m and  the  temperature  of  fusing  is, 
Tpusn  = 74  °C.  At  an  instant  of  time  during  a fire  scenario,  make  the  following  assumptions:  Tj^ 
= T,  = 70  °C;  u = 4.0  m/s  and  TG  = 1 10  °C  near  the  link;  throughout  the  space  of  fire  origin  the  gas 
is  still  relatively  transparent  and  does  not  participate  in  any  significant  radiative  heat  transfer 
exchange;  the  link  exchanges  radiation  with  far  field  surfaces  which  are  still  at  near-ambient 
temperature,  T^g  = 20  °C;  effective  emissivity  of  the  link  and  all  far-field  surfaces,  e,  are  0.9.  Then, 
using  Eq.  (7),  the  instantaneous  components  of  convective  and  radiative  heat  transfer  to  the  link  are 

Qconv.i  = h,(TG  - Tj)  = 1.72  x 103  W/m2 

4RAD.1  = a(T.tMB  - T})/(2/€  - 1)  = - 0.30  x 103  W/m2 
where  a = 5.67  x 10'8  W/(m2K4)  is  the  Stefan  Boltzman  constant. 

In  the  above  example,  it  is  seen  that  the  magnitude  of  radiation  cooling  of  the  link  to  the  far  field  is 
relatively  small  compared  to  that  of  convective  heating,  and  the  traditional  neglect  of  radiation  in  an 
engineering  calculation  of  link  response  may  be  justified.  However,  if  a higher-Tp^  link  is  used, 
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and  if,  say,  TLINK  = Tj  = 90  °C  and  all  of  the  other  above  parameters  are  the  same,  then 

qC0Nv,i  = 0-86  x 103  W/m2;  qRAD , = - 0.46  x 103  W/m2 

In  the  latter  case,  an  analysis  that  neglects  radiation  compared  to  convection  would  be  unacceptable. 
Also,  in  the  case  of  the  physically  larger  bar,  the  third  element  of  the  three-element  link-mount, 
Lconv.3  ~ 0.076  m > LC0NV ,.  Therefore,  From  Eq.  (7)  it  is  clear  that  h3  will  be  significantly  smaller 
than  hj,  and  it  can  be  anticipated  that,  compared  to  convection,  radiative  heating  of  this  bar  element 
will  always  be  important. 

Finally,  in  simulating  the  thermal  response  of  each  of  the  three  link-mount  elements  during  a fire 
scenario,  direct  radiation  from  the  combustion  zone  will  lead  to  an  additional  component  of  heat 
transfer.  This  can  lead  to  significant  modifications  in  the  response  of  all  the  elements,  i.e.,  even  for 
the  physically-small  link.  For  example,  for  a fire  with  total  energy  release  rate  QHRE  = 10  MW,  if 
35  percent  of  the  energy  released  in  the  combustion  zone  is  radiated  away  as  from  as  a point  source 
of  radiation,  then  the  rate  of  radiation  flux  incident  on  the  projection  of  a target  area  (e.g.,  the  surface 
of  the  fusible  link)  normal  to  a radius  from  the  source,  where  the  target  is  R = 10  m from  the  fire, 
will  be 


qRADFiRE  = 0.35Qhre/(47tR2)  = 2.78  x 103  W/m2 

To  estimate  the  significance  of  this  qRADFiRE  relative  to  qb0Nv,i  m the  present  example,  assume  a 
ceiling  height  of  7.5  m,  i.e.,  the  link  is  near  the  ceiling,  at  a radius  of  6.6  m from  the  point  of 
plume/ceiling  impingement.  Then,  at  the  radial  location  of  the  link  and  for  an  expansive  ceiling  with 
a thin  upper  layer,  the  fire  model  LAVENT  [6] [13] [14]  was  used  to  estimate  the  maximum  ceiling 
jet  velocity  and  temperature,  which  were  found  to  be  u = 2.1  m/s  and  TG  = 246  °C,  respectively. 
Assuming  again  that  TLINK  = Tj  = 70  °C,  Eq.  (7)  now  leads  to  the  result 

qC0Nv,i  =h1(Tc-T1)  = 2.10x  103  W/m2 

Thus,  in  this  example,  it  is  concluded  that  direct  radiative  heat  transfer  from  the  fire  to  the  link,  being 
comparable  to  convective  heat  transfer,  must  be  taken  into  account  in  an  analysis  of  the  thermal 
response  of  the  link.  It  therefore  seems  that,  in  general,  qRADFIRE  should  be  taken  into  account  in 
engineering  analyses  of  fusible-link  response.  Indeed,  except  in  the  case  of  particularly  small- 
dimension  fusible  links  where  pre-actuation  exposure  to  direct  fire  radiation  is  not  particularly 
intense,  the  latter  radiation  component  should  also  be  accounted  for  when  modeling  the  thermal 
response  of  the  fusible  links  used  to  actuate  sprinklers. 

In  view  of  the  last  observation,  it  is  noteworthy  that  in  real  fire  scenarios,  and  for  the  purpose  of 
making  a priori  estimates  of  qRADEiRE’  it  is  typically  impossible  to  estimate  accurately  the  effective 
location  of  the  combustion  zone,  both  in  terms  of  assigning  radial  distance  from  and  orientation  to 
the  radiated  surface  of  a target  fusible  link  or  other  link-mount  elements.  Since  qRADF1RE  is  strongly 
dependent  on  the  latter  two  variables,  the  situation  for  accurate  a priori  simulation  estimates  of  link 
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thermal  response  is  particularly  problematic.  Also,  to  the  extent  that  (in  the  relatively  early  time 
intervals  of  interest)  the  ceiling  jet  and  the  developing  upper  gas  layer  are  not  relatively  transparent 
to  relevant  radiation  wave  lengths,  the  above  radiative  heat  transfer  estimates  and  the  model 
equations  to  be  developed  below  would  have  to  be  modified  appropriately. 

In  order  to  include  the  above  effects  of  radiative  heat  transfer  on  element  response,  and  with  the 
latter  observation  in  mind,  model  equations  for  qRAD  s and  qRADnRE  wiU  be  constructed  consistent 
with  the  following  definition  and  assumptions: 

Prior  to  the  time  of  link  fusing,  assume  that  the  lower  surfaces  of  the  room  of  fire  origin  are 
still  near  ambient  temperature  and  that  radiation  absorption  and  emission  of  the  upper  gas 
layer  and  ceiling  jet  gases  are  negligible. 

Define  AC0NV  j as  the  effective  area  of  element  i to  which  gas-to-element  convective  heat 
transfer  takes  place. 

Regarding  qRADi,  assume  half  of  AC0NVi,  the  “top  " half  exchanges  radiation  with  afar 
field  having  temperature  and  effective  emissivity  identical  to  that  of  the  nearby  ceiling 
surface  of  the  room,  TCEIL  and  €CEIL,  respectively.  Assume  the  other  half  of  AC0NV 
the  “ bottom ” half,  exchanges  radiation  with  afar  field  having  temperature  and  effective 
emissivity  identical  to  that  of  the  floor  surface  of  the  room  or  of  the  “facing  ” surfaces  of  the 
room  contents,  T^g  and  €floor>  respectively,  i.e.,for  times  of  interest,  the  floor  or  facing 
content’s  surfaces  are  assumed  to  still  be  close  to  ambient  temperature. 

Regarding  q^on^  assume  that  the  radiant  energy  from  the  fire ’s  combustion  zone, 
effectively  a point  source  at  elevation  yFIRE  above  the  floor,  radiates  isotropically  at  a rate 
^radQftre’  where  Qhfie  is  the  total  energy  release  rate  of  the  fire  and  A RAD  is  the  fraction 
of  this  radiated  to  the  far  field,  of  the  order  of  0.35  for  typical  hazardous  fires  [6]. 

Assume  half  of  AC0NVj  is  exposed  to  direct  radiation  from  the  fire,  i.e.,  one  side  of  the 
element  is  illuminated  by  the  combustion  zone  and  the  other  side  is  shaded.  Assume  further, 
that  the  projection  of  this  exposed  surface  on  a surface  normal  to  a radial  vector  from  the 
fire  source  is  half  of  its  full  area. 

Model  equations  for  the  thermal  response  of  the  fusible  link.  From  the  onset  of  a fire  scenario, 
t = 0,  when  all  T,  are  at  T^g,  to  the  time  of  link  fusing,  t = tp^g  when  Tj  = TpuSE,  and  consistent 
with  the  above  definitions  and  assumptions  as  depicted  in  Figure  4,  energy  conservation  for  the  three 
elements  of  the  link-mount  leads  to 
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For  0 < t < tFUSE  = specified: 


mjCjdT/dt  = h,A( 


l^CONV.l 


m2c2dT2/dt  = h2A( 


^CONV^' 


■2(Tg  - T2)  + C - T2)  - C 23(^2  “ T3)  + QraD,2  + 0] 


RADHRE.2 


m3c3dT3/dt  = h 


(8) 


QraD.i  — (■^CONV,i^)^[(r^CEIL  * T|)/(1/€cEIL  + ~ 1) 


+ (T^B-Tt)/(l/eFLOoR+l/e,-l)]  (9) 


QraDFIRE.i  ~ (^CONV,/4)AradQfire  /{ 4tu [r  + (yCEIL  - Yhre)  ] } 


where 


T,(t  = 0)  = Tamb,  i = 1,2,  or  3 


(10) 


Qfire  — Qfire(I)’  Tg  — TG(t),  u — u(t),  and  TCEIL  — TCEIL(t)  are  specified 


In  Eqs.  (8),  the  hjA^-Qj^  j(TG  - Tj)  terms  represent  the  rate  of  convective  heat  transfer  from  the  gas  to 
element  i and  the  C^fT;  - Tj)  terms  represent  the  rate  of  conductive  heat  transfer  from  element  i to 
element  j.  In  Eqs.  (8)  and  (9),  the  Qrad  s and  Qradfire.i  terms  are  the  respective  total  rates  of  heat 
transfer  to  element  i from  radiation  exchanges  with  the  bounding  surface  of  the  enclosure  and  from 
direct  radiation  from  the  fire;  el  is  the  effective  emittance/absorptance  of  the  assumed-grey  surfaces 
of  element  i;  yCEIL  is  the  elevation  of  the  ceiling  above  the  floor;  and  r is  the  radial  distance  of  the 
link-mount  assembly  from  the  point  of  fire  plume/ceiling  impingement. 

Following  nomenclature  that  corresponds  to  and  extends  that  of  References  [15]  and  [16],  define 


tj  = miC/thjAcoNv  -),  i = 1,  2,  or  3 


i^CONV,i 


(ID 


and  rewrite  Eqs.  (8)  as 


dT/dt  — (Tg  - - C i2(Tj  - T2)/(m,Cj)  + FRAD  j 


(81 


dT3/dt  — (Tg  - T3)/t3  + C 23(^2 _ T3)/(m3c3)  + FRAD3 


where  the  F^  i?  which  include  all  effects  of  radiation  heat  transfer  on  element  i,  are 
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pRAD.i  — pRAD,i  (Tt;  specified  parameters)  — (Qradj  + QRADFiRE.iV(mici)  (12) 

Using  Eq.  (7)  in  Eq.  (11)  leads  to 

t,u1/2  - RTIi  - [(LC0NV j 1/2m,Ci/ AC0NV  ;)/2 . 96]  [m2K/(Ws1/2)]  (13) 

As  can  be  seen,  the  right  hand  side  of  the  above  equation  is  a constant  determined  by 
mechanical/thermal  and  geometric  properties  of  element  i.  For  the  fusible  link,  element  1,  this 
constant  is  identical  to  the  RTI  of  the  link  as  defined  in  Reference  [15].  Here,  the  right  hand  side 
of  Eq.  (13)  is  designated  as  RTIj,  or  the  RTI  of  element  i. 

Using  Eq.  (13)  in  Eqs.  (8')  leads  to 

For  0 < t < tFUSE  = specified: 

dT/dt  = [T0  - (1  + C12/u1/2)Ti  + (C12/u1,2)T,]ui,2/RTIi  + FRAD , (8") 

dT2/dt  = [TG  - (1  + C2|/u1'2  + Cnlam)l2  + (C2I/u1,2)T,  + (C23/ul,2)T3]u1,2/RTI2  + F^, 
dTj/dt  = [T0  - (1  + C32/u,,2)T3  + (C32/ui,2)T2]u1/2/RTI3  + FRAD3 
subject  to  the  conditions  ofEq.  ( 10) 


where 


C12  - C'j.RTIj^mjCj)  = C '12[(LC0NV  j 1,2 1 AC0NV , )/2.96] [m2K/(W s1/2)] ; 
C21  ^ C"12RTI2/(m2c2)  = C '12[(LC0NV 21/2/ AC0NV 2)/2.96] [m2K/(W s 1/2)] 
C23  - C"23RTI2/(m2c2)  = C'23[(LCONV21/2/ACONV2)/2.96][m2K/(Ws1/2)]; 
C32  - C"23RTI3/(m3c3)  = C '23[(LCONV  3 1/2/ AC0NV  3)/2 .96]  [m2K/(Ws 1/2)] 
Also,  from  Eqs.  (14) 


C 21=  C12(m,c1RTI2)/(m2c2RTI1); 

(15) 

C32=  C23(m2c2RTI3)/(m3c3RTI2) 

Note  that  m,  and  c,  can  be  easily  obtained.  Therefore,  it  is  assumed  they  are  known  and  specified. 
Then,  corresponding  to  a set  of  five  system  parameters  RTIL  RTI2.  RTI3,  C12,  and  C23,  to 
be  determined,  the  values  of  C21  and  C32  are  obtained  directly  from  Eqs.  (15). 
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Determining  the  Thermal  Properties  of  the  Link-Mount  Assembly 


General  considerations.  The  solution  to  the  problem  of  Eqs.  (8")  and  (10)  would  provide  a 
prediction  of  the  thermal  response  of  the  link-mount  assembly  when  it  is  exposed  to  a specified  fire 
environment  and  corresponding  local  flow.  For  a vent  located  at  a particular  location  relative  to  a 
fire  scenario  of  interest,  and  with  known  local  values  of  TG(t)  and  u(t)  (e.g.,  determined  from  a fire 
model  simulation),  the  solution  would  yield  the  Tj(t).  Of  interest  here  would  be  the  solution  for  T (t), 
the  temperature  of  the  fusible  link  itself,  and,  in  particular,  a determination  of  tp^,  the  time  when 
it  fuses  and  when  the  vent  opens. 

To  fully  define  the  problem  for  a particular  vent  design,  it  is  necessary  to  determine  the  thermal 
properties  of  the  link-mount  assembly,  i.e.,  the  five  system  parameters  RTIL  RTI2.  RTI3,  C12,  and  C23. 
This  section  addresses  means  of  determining  these  parameters.  This  will  use  a combined- 
experimental/analytic  method,  where  the  experimental  aspect  involves  use  of  the  plunge  test  [1],  [2], 
[15],  originally  developed  to  determine  the  RTI  for  sprinkler  links. 

The  plunge  test.  The  plunge  test  is  a wind-tunnel  test  where,  in  its  designed  application,  a sprinkler 
head,  with  its  actuating  link-mount  assembly  (Tp^  of  the  link  is  known),  all  at  initially  ambient 
temperature,  is  inserted  (i.e.,  “plunged”)  into  an  air  stream  of  specified  constant  temperature,  TG  > 
Trise’  and  specified  constant  flow  speed,  u.  An  original  purpose  of  the  plunge  test  was  to  determine 
the  thermal  characteristics  of  different  sprinkler  designs,  so  that  sprinkler  actuation,  i.e.,  fusing  of 
the  link,  could  be  predicted  when  the  sprinkler  head  is  exposed  to  an  arbitrary,  specified,  time- 
dependent  fire-environment. 

Note  that  the  currently-accepted  model  equation  for  sprinkler  link  thermal  response  (see,  e.g., 
References  [15]  and  [16]),  which  includes  the  effect  of  conduction  loss  from  the  link  to  the  sprinkler 
base,  is  a special  case  of  Eqs.  (8")  and  (10),  where,  as  in  the  present  analysis,  T]  is  the  temperature 
of  the  link,  but  where  T2  is  the  temperature  of  the  base  of  the  sprinkler  assembly,  assumed  to  be  at 
a known  constant  temperature,  usually  taken  to  be  T^g;  C12  characterizes  the  conduction  path 
between  the  link  and  the  sprinkler  base;  radiation  effects  are  ignored,  and  only  the  first  of  Eqs.  (8") 
is  applicable.  Thus,  the  traditional  model  equations  for  the  response  of  a sprinkler  link  exposed  to 
an  arbitrary  fire  environment  are 

For  0 < t < tpusg  = specified: 

dT/dt  = [Tc  - (1  + C|2/u1,2)T,  + (Qj/OTJu^/RTI,  (8"') 

where 

Tj(t  = 0)  = T^g  and  TG  = TG(t),  u = u(t)  are  specified  (10') 

In  the  case  of  a constant-velocity/temperature  plunge-test  exposure,  and  neglecting  effects  of 
radiation,  the  general  solution  to  the  above  for  Tj(t)  is  readily  determined  (see  Eq.  (21)  of  Reference 
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[15]).  In  particular,  at  t = tpusE 


RTIj  — - tFUSEU  (1  + ^-'12^  " (TfuSE  “ TAMb)0  + ^12^U  VO-G  ^AMb)]  0^) 

Thus,  if  C12  is  known,  the  RT^  of  the  sprinkler  link  can  be  determined  from  a single  plunge  test. 


As  discussed  in  Reference  [15],  C12  can  be  determined  from  additional  application  of  the  plunge  test 
apparatus.  The  idea  is  to  expose  the  sprinkler  head  to  series  of  quasi-steady  exposures  at  a fixed  TG. 
The  exposures  involve  an  initial  value  of  u that  is  less  than  required  to  affect  fusing  of  the  link.  A 
steady-state  condition  is  established  where  the  conductive  heat  transfer  from  the  link  to  the  cool 
sprinkler  base  is  large  enough,  and  the  rate  of  heat  transfer  from  the  relatively  low-speed  gas  to  the 
link  is  small  enough  to  maintain  the  link  at  a temperature  less  than  Tp^.  Then  u is  increased  in 
small  steps.  After  allowing  the  system  enough  time  between  steps  to  reach  a quasi-steady  state,  with 
dT/dt  effectively  zero  (in  the  tests  of  Reference  [15],  an  interval  of  10  min  was  adequate),  the  value 
of  Tbase,  the  temperature  of  the  base  of  the  sprinkler,  is  recorded.  Eventually,  during  one  of  these 
steps,  the  link  will  attain  Tpusg.  The  u for  this  step  is  designated  as  the  critical  u,  uCR  = uCR(TG). 
Then,  the  steady- version  of  the  first  of  Eqs.  (8")  leads  to  the  desired  objective 


Ci2  — UCR  (Tc  - TpusE 


CR 


)/(TpusE  Tbase) 


(17) 


(also  see  Eq.  (25)  of  Reference  [15]),  where,  in  these  sustained  quasi-steady  exposures,  TBASEmay 
be  measurably  higher  than  T^g. 


Having  determined  C12  experimentally  from  the  latter  test  sequence,  it  can  be  used  in  Eq.  (16), 
together  with  the  previously  measured  plunge-test  value  of  tp^,  to  determine  the  desired  value  of 
RTIj.  With  the  values  of  RTIj  and  C12  in  hand,  the  response  of  the  link  to  any  arbitrary,  specified 
fire  exposure  can  be  determined  from  a solution  to  Eqs.  (S'")  and  (10"). 

The  plunge  test  and  the  link-mount  assembly.  As  is  the  case  for  single  sprinkler-head  assemblies, 
the  response  to  plunge  test  exposures  of  the  more  complex  link-mount  assembly  of  current  interest 
can  also  be  measured.  However,  while  the  above  two-stage  test  approach  (a  single  time-dependent 
exposure  to  link  fusing,  and  a series  of  quasi-steady  exposures)  leads  to  a determination  of  RT^  and 
C12,  and  a complete  thermal  characterization  of  the  two-parameter  sprinkler-head  systems,  this  is  not 
the  case  for  the  five-parameter  link-mount  assembly.  Futhermore,  because  of  the  current  concern 
that  radiation  effects  may  be  important,  even  in  the  plunge  test  environment  (i.e.,  radiation  exchanges 
between  the  elements  and  the  relatively  high-temperature  tunnel  walls),  the  situation  is  complicated 
even  farther.  In  short,  an  alternative  method  is  required  to  determine,  from  measured  time- 
dependent  responses  to  plunge  test  exposures,  the  five  system  parameters  of  the  link-mount 
assembly. 


Modified  Equations  for  Response  Under  Plunge  Test  Conditions 

The  problem  of  Eqs.  (8"")  and  a simplified  version  of  Eqs.  (10)  can  be  used  to  simulate  the  response 
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of  the  link-mount  assembly  during  a plunge  test  exposure  if  the  radiation  terms,  i?  as  defined  by 
Eqs.  (9)  and  (12),  are  modified/reinterpreted.  First,  the  non-operative  direct-fire-radiation  term, 
Qradfire.p  is  eliminated.  Then,  in  Q^AD , , the  characteristics  of  the  ceiling  and  the  floor  of  the  fire 
room  are  replaced  by  the  corresponding  characteristics  of  the  inside  wall  in  the  test  section  of  the 
tunnel,  which  surrounds  the  assembly.  This  leads  to 

FraD.i  = FRAD  i PLUNGE  = (QRAD,iV(miCi)  = AcONV,i°(TWALL  " Tj)/[(  INWALL  + l^i  _ 1 XmiCi)] 

(18) 

This  is  used  in  Eqs.  (8"),  and  the  problem  definition  is  completed  with 

Tj(t  = 0)  = T^g,  i = 1,  2,  or  3;  TG,  u,  and  TWALL  are  specified  constants  (10"') 


Determining  the  Thermal  Properties  of  the  Link-Mount  Assembly 

An  algorithm  for  determining  the  system  parameters.  Corresponding  to  the  unknown  RTI/s  and 
the  Cjj’s,  of  the  link-mount  assembly,  define  the  five  system  parameters  A,  by 

A,  = Tgu1/2/RTIj;  A,  = Tgu1/2/RTI2;  A3  = Tgu1/2/RTI3; 

(19) 

A4  = (C12/RTIj);  A5  = (C23/RTI2) 

For  an  arbitrary  choice  of  the  parameters,  let  T,(t;  ...,  A5)  be  the  solution  to  Eqs.  (8~),  (18)  and 

(10"),  from  t = 0 to  t = tp^.  Such  a solution  can  presumedly  be  obtained  by  numerical  integration. 

Let  the  link-mount  assembly  be  subjected  to  a plunge  test  and  denote  the  experimental  data  for  the 
variables  T,(t)  as  TEXPii(tj),  j = 1 , NF,  where  NF  is  the  number  of  times  that  data  were  acquired  up  to 

^FUSE- 

Now  consider  the  following  problem 

Determine  the  values  of  Ak  that  minimize 

3 Nf  (20) 

G(Ak,  TEXPi(tj))  = E E [TEXPi(tj)  - T j ( t j ; A1? ...,  A5)]2 

i=i  j=i 

It  is  now  assumed  that  the  set  of  Ak  solutions  to  the  problem  of  Eq.  (20)  correspond  to  a best  estimate 
of  the  desired  thermal  properties  of  the  link-mount  assembly.  Thus,  from  Eqs.  (19)  the  desired  set 
of  system  parameters  are  obtained  from 

RTIi  = TqU^/Aj,  i = 1,  3;  C12  = TGul,2(A4/A,);  C23  = TguI/2A5/A2  (21) 
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Characterization  of  the  system  is  completed  by  using  Eq.  (15)  and  the  results  of  Eq.  (21)  to  determine 

C21  = Tgu1/2(A4/A2)(m1c1)/(m2c2);  C32  = Tgu1/2(A5/A3)(m2c2)/(m3C3)  (22) 

Implementing  the  least  squares  curve  fitting  algorithm.  An  implementation  of  the  above 
algorithm  requires:  plunge  test  data,  TEXPi(tj);  numerical  integration  of  the  initial  value  problem  of 
Eqs.  (8"),  (18),  and  (10"),  necessary  to  compute  Tj(t;  Ak)  and  G(Ak,  TEXPi(tj))  of  Eq.  (20);  and  an 
algorithm  and  associated  computer  software  to  find  the  set  of  parameters  that  minimize  the  sum  of 
squares  of  a set  of  non-linear  functions  of  several  variables. 

In  the  carrying  out  the  above,  the  numerical  integration  subroutine  used  was  ODEINT  from 
Reference  [17],  and  the  minimization  subroutine  used  was  SNLS1E  from  Reference  [18]. 

Plunge  Tests  on  the  Link-Mount  Assembly  of  Figures  2 and  3 and  on  Related  Link-Mount 
Assemblies 

A modified  short-bar  link-mount  test  assembly.  As  mentioned,  the  bar  element  of  the  link-mount 
assembly  is  1.2  m long,  too  long  to  be  tested  in  the  0.20  m x 0.20  m test  section  of  the  plunge  test 
apparatus.  However,  it  is  expected  that  somewhat  away  from  its  end  mounts,  the  thermal  response 
of  this  full-length  bar  to  exposures  of  Figure-2-type  fire  environments  does  not  vary  significantly 
along  its  length.  It  is  also  reasonable  to  expect  that  the  thermal  response  to  fire  environments  of  a 
relatively-short  “free-end”  segment  of  the  bar  (length  of  at  least  a few  times  the  characteristic 
dimension  of  the  bar’s  section)  would  be  substantially  similar  to  that  of  the  full-length  bar. 

The  above  expectations  are  the  basis  for  testing  a modified  Figure  2-type  three-element  link 
assembly,  with  a shortened  free-end  bar  element,  and  for  assuming  that  its  response  in  both  uniform, 
steady,  plunge-type  flows  and  in  time-dependent  real  fire  environments  will  be  substantially  similar 
to  that  of  a full-length-bar  assembly.  The  length  of  the  bar  element  of  the  tested  assembly  was  0.121 
m.  Except  for  the  shortened  length  of  the  bar,  the  construction  of  the  three-element  link-mount  test 
assembly,  i.e.,  the  tab  and  link  elements  and  the  means  of  element-to-element  attachment,  per  Figure 
3,  was  identical  to  that  found  in  the  real  vent  hardware. 

The  fusible  links  used  were  of  all  of  identical  design,  as  depicted  in  Figures  1 and  3,  with  a rated  fuse 
temperature  of  74  °C. 

The  test  procedure  and  measurements.  See  Figure  5.  The  test  assembly  was  mounted  via  an 
insulating  0.010  m-diameter  stand-off  rod  between  the  base  of  its  bar  element  and  a steel  cover  plate 
used  to  cover  the  opening  of  the  test  section.  The  design  led  to  rigid  test  fixture  with  a 0.020  m 
stand-off  distance  between  the  bar  base  and  the  cover. 

In  the  test  procedure,  the  tunnel,  with  covered  ceiling  test  section  opening,  but  without  the  test 
fixture,  was  brought  up  to  the  desired  steady  gas  velocity  and  temperature;  the  inside  wall 
temperature,  continuously  measured  in  the  vicinity  of  the  test  section,  was  also  at  steady  state.  At 
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the  initiation  of  a test,  the  test  section  was  opened,  and  the  test  assembly  fixture  was  inserted  upside- 
down  into  the  test  section.  The  insertion  action  took  on  the  order  of  1 s. 

The  temperature  of  each  of  the  three  elements  was  measured  every  0.5  s with  thermocouples  peened 
to  the  element  surface  at  the  locations  indicated  in  Figure  5-  twall  was  measured  in  the  immediate 
vicinity  of  the  test  section  and  TG  was  measured  upstream  of  the  test  section. 

One  test  was  carried  out  with  the  air  flow  normal  to  the  length  of  the  bar  axis  (as  in  Figure  5)  and 
one  test  with  the  air  flow  parallel  to  the  bar  axis. 

A feature  of  the  assembly/cover  test  fixture,  indicated  in  Figure  5,  allowed  for  a constant  1.4  Kg 
axial  spring  loading  of  the  link,  and  for  visual  confirmation,  outside  the  tunnel,  of  when  the  link 
fused.  The  time  of  link  fusing  could  also  be  clearly  distinguished  by  a significant  subsequent 
increase  of  temperature  of  the  link  thermocouple  which  quickly  approached  the  air  flow  temperature. 

A test  was  terminated  at  t = tp^. 

Related  link-mount  assembly:  three-element  assembly  with  insulating  connections.  One  test 

was  carried  out  on  a second  three-element  assembly,  referred  to  here  as  the  “element-isolated 
assembly.”  This  was  identical  to  the  three-element  assembly  of  Figure  3,  except  that  it  was 
fabricated  so  that  the  element-to-element  thermal  conduction  paths  are  modified  in  such  a way  as 
to  significantly  reduce  inter-element  conduction,  without  modifying  the  basic  gas-to-element 
convective  heat  transfer  phenomena.  In  the  assembly  design,  the  link-to-tab  and  tab-to-bar  joints 
involved  insulating  spacers  fabricated  from  0.0016  m-thick  epoxy-resin-type  circuit  board  material 
and  fastened  to  the  respective  surfaces  with  a high-temperature  adhesive. 

The  test  on  the  element-isolated  assembly  was  with  the  air  flow  parallel  to  the  test  assembly. 

Related  link-mount  assembly:  simply  supported  fusible  link.  A final  test  was  carried  out  on  a 
simply-supported  fusible  link  (see  Figure  1).  As  in  the  normal-flow  test  of  Figure  5,  the  surface  of 
the  link  was  normal  to  the  flow,  and  it  was  instrumented  with  a thermocouple  peened  to  the  center 
of  its  rear  surface.  The  test  fixture  provided  spring  loading  of  the  link  during  the  test  and  visual 
indication  of  link  fusing. 

Test  conditions  and  test  results.  For  all  tests,  thermostatic  control  of  the  gas  flow  temperature  at 
135  °C  led  to  measured  values  of  TG  = 134  °C  ± 1.0  °C.  Using  previous  calibrations  of  tunnel  flow 
velocity  V5  gas  flow  temperature,  the  tunnel  fan  speed  control  was  set  at  a value  corresponding  to  an 
average  test  flow  velocity  of  u = 2.5  m/s. 

Table  2 presents  an  overview  of  all  tests,  including:  test  conditions;  the  test  results:  TWALL, 
T,(t  = tpusp),  tpuSE,  Nf;  and  reference  to  the  figure  numbers  where  the  measured  T,  are  plotted.  Also 
included  in  the  table  are  tMELT  and  NM,  to  be  defined  below. 
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Table  2. 

Plunge  test  configurations,  conditions,  and  results. 

test 

no. 

configuration 

flow 

direction 

TaMB 

tg 

[°C± 

T\VALL 

1°C] 

T infuse) 

fnjsE 

[s] 

NF 

fviELT 

[s] 

nm 

Fig. 

no. 

1 

full  assembly 

normal 

25. 

134. 

90. 

74. 

68.0 

137 

34.1 

69 

6 

2 

full  assembly 

parallel 

24. 

134. 

97. 

74. 

63.0 

127 

34.5 

70 

7 

3 

element- 

isolated 

assembly 

parallel 

24. 

134. 

97. 

73. 

54.5 

110 

31.5 

64 

8 

4 

simply- 

supported 

link 

normal 

24. 

133. 

100. 

74. 

53.5 

108 

25.0 

50 

9 

Overview  of  Results 

Main  features  of  the  test  results  are: 

To  within  experimental  accuracy,  the  links  used  in  each  of  the  five  tests  fused  at  the  specified  device 
fuse  temperature,  T^tp^g),  of  74  °C. 

Results  of  tests  1 and  2 indicated  that  orientation  of  the  three-element  assemble  relative  to  flow 
direction  does  not  have  a significant  effect  on  tp^p  (68.0  s and  63.0  s,  respectively)  or  on  the  T,  (see 
Figures  6 and  7,  respectively). 

From  Figures  6-8,  the  effect  of  increasing  element-to-element  thermal  isolation  (compare  the  full 
assembly  of  tests  1 and  2 to  the  element-isolated  assembly  of  test  3)  is  seen  to  have  qualitatively 
modified  element  response.  Initially,  the  temperature  of  both  the  bar,  T3,  and  the  tab,  T2,  in  test  3 
increased  more  rapidly  than  they  had  in  tests  1 and  2,  with  T2  greater  than  the  temperature  of  the  link, 
Tj,  until  well  into  the  test. 

It  is  seen  that  the  differences  between  tp^  for  the  full  assembly  in  tests  1 and  2 (68.0  s and  63.0  s, 
respectively)  and  tpuSE  for  the  simply  supported  link  of  test  4 (53.5  s)  were  14.5  s and  9.5  s, 
respectively,  a relative  reduction  of  21  % and  15  %,  respectively.  Thus,  in  spite  of  the  firm 
connection  of  the  link  to  the  relatively  massive  tab  and  bar  in  the  three-element  link  mount  assembly, 
tpusE  for  the  full  assembly  is  not  increased  significantly  over  that  of  the  simply  supported  link. 
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In  terms  of  a tpuSE  measure,  the  response  of  the  link  of  the  element-isolated  assembly  (54.5  s)  in  test 
3 was  very  close  to  that  of  the  simply  supported  link  (53.5  s)  of  test  4. 

For  all  tests,  the  Tj  histories  of  Figure  6-9  indicate  that  the  time  interval  over  which  link-solder 
melting  takes  place,  i.e.,  the  interval  when  the  link  is  maintained  close  to  its  fuse  temperature,  is  a 
significant  fraction  of  tFUSE,  of  the  order  of  fifty  percent. 

Relative  to  the  success  of  the  basic  RTI  modeling  concept,  which  neglects  the  time  interval  of  solder 
melting  compared  to  tp^,  the  latter  result  is  particularly  significant.  Thus,  for  the  fusible  link  design 
tested  here,  the  RTI  model  equation  for  the  link,  whether  for  a simply  supported  link,  i.e.,  Eq.  (8"'), 
or  for  the  three-element  link  mount  of  current  interest,  i.e.,  the  first  of  Eqs.  (S'"),  or  of  Eqs.  (8"), 
etc.,  is  not  adequate.  Rather,  the  model  equations  will  have  to  be  modified  to  take  into  account  the 
phenomenon  of  solder  melting  and  the  extended  time  of  nearly-constant  link  temperature  prior  to 
link  fusing. 

Applying  the  Least  Squares  Curve  Fitting  Algorithm  to  Determine  System  Parameters 

Two  applications  of  the  algorithm  for  each  test.  With  the  availability  of  the  above  T,  test  data, 
it  is  now  possible  to  apply  the  above-proposed  least  squares  algorithm  to  determine  the  RTI;  and  the 
CI:I  parameters  that  define  the  thermal  response  characteristics  of  the  three  configurations:  the  full 
assembly  (tests  1 and  2),  the  element  isolated  assembly  (test  3)  and  the  simply  supported  link  (test 
4).  This  will  be  done  without  regard  for  the  above  observation  on  the  inadequacy  of  the  model 
equations  for  Tj.  Thus,  although  it  is  now  apparent  that  the  previously  adopted  system  of  equations 
does  not  represent  an  optimum  model  for  the  Ti?  its  is  still  reasonable  to  use  the  experimental  plunge 
test  results  and  to  determine  those  values  of  RTI,  and  Qj  that  optimize  the  system  equations,  Eqs. 
(8")  and  (10"),  according  to  the  algorithm  of  Eq.  (20),  etc.  Later,  overall  model  equations  will  be 
improved,  and  corresponding  system  parameters  will  be  determined  with  a modified  methodology 
that  takes  into  consideration  the  above  link-melting  phenomenon. 

The  algorithm  was  applied  twice  for  each  of  the  four  tests.  The  first  application  involved 
experimental  data  from  the  initiation  of  the  test  up  to  the  time,  t = tMELT,  when  the  measured  link 
temperature,  TL  reaches  the  plateau,  close  to  Tfuse  = 74  °C,  that  clearly  indicates  initiation  of  the 
link-melt  process.  Note  that  in  this  0 < t < tMELT  time  interval,  the  original  equation  set  is  still 
expected  to  represent  a valid  model,  since  the  objections  raised  earlier  regarding  phenomenon  of 
solder  melting  has  not  yet  been  initiated.  As  can  be  seen  in  Figures  6-9,  there  is  no  ambiguity  in 
choosing  the  experimental  time  and  corresponding  T,  sets  that  correspond  to  t = tMELT.  For  each  test, 
define  NM  as  the  number  of  times  that  data  were  acquired  up  to  tMELT.  Values  of  tMELT  and  NM 
are  included  in  Table  2. 

For  each  of  the  tests,  the  second  application  of  the  algorithm  involved  all  of  test  data  in  the  interval 
0 < t < tEUSE. 

An  initial  estimate  of  the  system  parameters.  The  essential  subroutine  of  the  algorithm  software 
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is  SNLS1E.  To  use  this,  it  is  necessary  to  provide  as  input  an  initial  estimate  of  the  system 
parameters.  Based  on  Eqs.  (8),  the  geometry  of  the  elements,  their  material  properties, 
measurements  of  their  masses  (nq  = 0.0156  kg;  m2  = 0.0280  kg;  m,  = 0.780  kg),  and  using  Eqs. 
(13)  and  (14),  an  initial  estimate  for  the  parameters  of  the  three-element  link  mount  assembly  was 
determined  to  be 

estimate  of  three-element  link  mount  parameters:  (23) 

RTIj  = 200  (msf;  RTI2  = 430  (msf;  RTI3  = 777  (ms)1/2;  C12  = 1.04  (ms)*;  C23  = 7.7  (ms)* 

The  optimum  “solution”  parameters  and  corresponding  simulated  responses.  Results  of  the 
above  algorithm  applications,  i.e.,  the  “best-fit”  sets  of  values  for  the  RTIj  and  C1J?  are  presented  in 
Table  3.  As  indicated  in  the  table,  simulations  for  the  test  Tj(t)’s  that  correspond  to  these  sets  of 
values,  i.e.,  the  solutions  to  Eqs.  (8"),  (18)  and  (10"),  are  presented  in  Figures  10  - 13  for  the  0 < 
t < tMELT  data  fits  and  in  Figures  14  - 17  for  the  0 < t < ^jSE  data  fits.  Included  in  Figures  10-17 
are  the  applicable  experimental  data  as  were  presented  previously  in  Figures  6-9. 

As  can  be  seen  in  Figures  10  - 13,  in  the  range  0 < t < tMELT  and  for  the  0 < t < tMELT  data  fits  there 
is  generally  excellent  agreement  between  the  experimental  data  and  the  model  simulations.  Thus, 
for  exposure  to  the  particular  plunge  test  environment  of  tests  1 - 4,  and  at  least  up  to  the  initiation 
of  link  solder  melting,  the  algorithm  was  successful  in  identifying  sets  of  system  parameters,  the  ones 
solved  for  and  presented  in  Table  3 under  the  heading  “fit  from  t = 0 to  t = tMELT,”  that,  when  used 
with  the  model  equation  set  Eqs.  (8"),  (18)  and  (10""),  simulate  well  the  measured  response  of  the 
system  configurations  tested.  Further,  it  is  reasonable  to  expect  that  the  same  model  equation  set, 
with  the  same  solution  parameters  of  Table  3,  but  with  arbitrary  Eq.  (10"")-conditions,  will  provide 
similarly  accurate  simulations  for  the  system  responses  to  arbitrary  plunge  test  environments. 
Additional  testing  is  required  to  determine  if  the  latter  expectation  is  realized.  In  response  to  the 
main  objective  of  this  work,  it  is  finally  expected  that,  when  used  with  the  general  equation  set  Eqs. 
(8""),  (9),  (10),  and  (12),  the  solution  system  parameters  of  Table  3 will  lead  to  accurate  simulations 
of  the  thermal  response,  in  the  range  0 < t < tMELT  , of  the  configurations  to  arbitrary  Eq.  (10)- 
specified  fire  environments. 

As  can  be  seen  in  Figures  10  - 13,  in  the  range  0 < t < tp^,  and  in  terms  of  predicting  the  critical 
event  of  system  response,  i.e.,  the  time  t = tFUSE  when  = 74  °C,  the  agreement  between  the 
experimental  data  and  the  model  simulations  is  not  good.  In  particular,  in  the  figures  and  for  the  four 
tests,  the  predicted  tpu^’s  are  found  to  be  relatively  25  % to  30  % less  than  the  corresponding 
measured  tpj^’s.  Thus,  as  suggested  earlier,  for  the  fusible  link  design  used  in  tests  1 - 4.  and,  no 
doubt,  for  other  relatively-large-solder-mass  link  designs,  the  governing  equation  for  Tj  has  to  be 
modified  so  that  it  models  adequately  the  near-constant-^  solder-melt  phase  of  link  response. 

Modifying  the  Model  Equations  to  Account  for  the  Solder-Melting  Phase  of  Link  Response 

Simulating  the  solder-melting  phase  of  link  response.  For  a link  design  with  arbitrary  specified 
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Table  3.  System  parameters  determined  from  the  least  squares  curve 
fitting  algorithm  (all  e’s  of  Eq.  (18)  assumed  to  be  0.9). 


fit  from  t = 0 to  t = tMELX  fit  from  t = 0 to  t = tFUSE 

test  number:  test  number: 


1 

2 

3 

4 

1 

2 

3 

4 

RTIjAmsf  = 

83.01 

76.53 

88.16 

71.96 

66.58 

65.54 

102.5 

52.17 

RTytmsf  = 

100.1 

46.40 

15.45 

NA 

119.7 

70.13 

19.27 

NA 

RTIj/Cms)'7- 

1488 

1907 

1907 

NA 

1448 

1547 

2293 

NA 

C12/(ms)1/2  = 

1.401 

1.684 

0.8286 

NA 

3.334 

3.060 

6.370 

NA 

Cyims)'72  = 

12.44 

14.04 

5.065 

NA 

13.83 

13.41 

4.409 

NA 

Total 

Residuals/N  = 

0.0261 

0.0296 

0.0767 

0.0340 

0.0675 

0.0584 

0.0962 

0.0880 

Figure  no. 

10 

11 

12 

13 

14 

15 

16 

17 

Tfuse,  the  near-constant-^  solder-melting  phase  of  link  response  can  be  closely  simulated  by  adding 
the  following  constraint  to  the  model  equations,  Eqs.  (8")  or  (S'"): 

if  Tj  > Tpusg,  then  dT/dt  = 0 (24) 

Note  from  the  first  of  the  model  equations,  the  ones  for  Tls  and  from  Eq.  (13),  where  it  is  seen  that 
RTI,  is  proportional  to  cjs  that  the  above  constraint  is  equivalent  to  an  assumed  arbitrarily-large  step 
increase  in  c1?  in  the  temperature  range  Tj  > Tp^,  from  its  specified  constant  low-temperature 
value.  This  increase  in  effective  specific  heat  of  the  link  simulates  the  action  of  the  latent  heat  of 
melting  of  the  solder,  which  maintains  the  entire  link  component  at  the  nearly-constant  temperature, 
Tfuse’  as  the  solder  absorbs  the  relatively  large  amount  of  energy  required  to  bring  it  from  a solid  to 
liquid  state. 

With  the  Eq.  (24)-modified  model  equations,  the  curve-fitting  algorithm  can  be  applied  again,  as 
above.  The  result  would  be  sets  of  solution  system  parameters,  to  be  used  to  revise  Table  3 entries 
under  the  heading  “fit  from  t = 0 to  t = tp^.” 
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It  is  clear  from  a visual  inspection  of  the  data,  e.g.,  Figures  14  - 17,  that  the  revised  model  equations, 
used  with  the  latter  revised  solution  system  parameters,  would  lead  to  improved  goodness  of  fit 
between  simulation  and  experiment,  in  the  sense  of  reducing  G of  Eq.  (20).  However,  with  regard 
to  providing  an  accurate  simulation  estimate  for  the  critical  value  tpuSE,  the  situation  is  far  worse  than 
before.  Thus,  while  the  original  simulation  estimates  for  tpuSE  were  25  % to  30  % less  then  the 
measured  values,  the  revised  model  equations,  i.e.,  with  the  constraint  of  Eq.  (24),  do  not  yield  any 
unique  prediction  of  tpuSE! 

To  resolve  the  situation,  additional  modeling  detail  of  the  solder-melting  process  is  required.  In  this 
regard,  a variety  of  different  modeling  approaches  are  possible,  e.g.,  an  approach  where  the  model 
accounts  for  time-dependent  changes  of  m]  due  to  loss  of  solder  (dripping)  during  the  melting 
process.  (In  general,  the  mass  of  the  pre-melt  solder  component  of  the  link  can  be  a significant 
fraction  of  m,.)  However,  to  achieve  the  stated  objective  of  this  work  it  seems  that  an  elaborate 
analysis  is  unnecessary. 

The  Seven-Step  Procedure  for  determining  system  parameters  and  establishing  final  model 

equations.  Consistent  with  the  above  results  and  ideas,  an  alternative  variation  to  the  Eq.  (24)- 
approach  will  now  be  presented  that  yields  a general  model  for  simulating  the  thermal  response, 
including  a prediction  for  tpuSE,  of  any  particular  three-element  link-mount  system  configuration. 

Values  for  the  system  parameters  are  determined  and  the  final  model  equations  are  established  by 
carrying  out  the  following  steps,  referred  to  below  as  the  “Seven-Step  Procedure:” 

1 . Acquire  plunge  test  data  for  the  system  configuration  to  be  modeled,  similar  to  the  data  of 
Tests  1 - 4,  where  the  configuration  includes  a particular  fusible  link  design  with 
manufacturer-specified  Tp^.  From  the  data,  identify  the  value  of  tpuSE. 

2.  Using  the  curve  fit  algorithm,  determine  system  parameters  RTI,  and  Cy  associated  with  the 
solution  to  the  model  equations,  Eqs.  (8~),  (18),  (10")?  and  (24),  that  yields  a best  fit  to  the 
data  in  the  range  0 < t < tp^.  In  the  latter  implementation  of  the  algorithm,  the  values  of  all 
e’s  of  Eq.  (18)  are  assumed  to  be  known  and  specified.  (For  results  reported  in  this  work, 
all  e’s  were  always  specified  as  0.9.) 

3.  Consistent  with  the  above  plunge  test  data,  select  a relatively-small  (compared  to  AT  = TpuSE 

- T^g)  temperature  increment,  ATmelt,  such  that  the  relatively-constant-temperature  solder- 
melt  region  of  the  Tj  data  of  the  test  of  item  1 are  in  the  range  TpuSE  - ATmelt  < < TpuSE. 

4.  Define  the  solid-to-liquid  phase  change  of  the  link  solder  as  being  initiated  at  temperature 
Tmelt’  where 

Tmelt  = TpusE  - ATmelt  (25) 

From  the  experimental  data,  identify  the  time  tMELT,  when  T,  first  reaches  Tmelt. 
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5.  Define  RTImelt  as  the  new,  effective,  melt-phase  value  of  RTIj,  to  be  determined 
below.  (This  will  have  the  property:  RTImelt  » RTIj.)  Corresponding  to  RTImelt 
is  cMELT,  the  newly-defined,  effective,  melt-phase  value  of  c, 

cmelt  “ CiRT^nLj/RT^  (26) 

Add  the  following  constraint  to  the  model  equations: 

if  tmelt  ^ Ti  < Tpusg,  then,  in  Eqs.  (8"),  (12),  and  (18), 

replace  RTl{  by  RTImelt  and  cl  by  cMELT  (27) 

6.  Find  the  value  of  RTImelt  that  leads  to  a solution  to  Eqs.  (8"),  (18),  (10"),  and  (27), 
integrated  in  the  range  0 < t < tpuSE,  that  satisfies  the  condition 

^1  (1-fuse)  = ^FUSE  (28) 

Designate  this  value  of  RTImelt  as  a new  fixed  system  parameter,  in  addition  to  the  RTI,  and 
Cy  determined  above  in  item  4. 

7.  The  final  model  equations  are:  Eqs.  (8"),  (9),  (12),  (10),  and  (27),  for  fire  scenarios;  or,  in 
the  case  of  tunnel  test  simulations:  Eqs.  (8"),  (18),  (10"),  and  (27).  In  general  applications 
of  the  model,  fusing  of  the  link  is  predicted  to  occur  at  the  first  simulated  time  when  Tj  = 

TrjSE- 

In  the  above  procedure,  step  2 insures  a good  general  fit  between  the  model  simulation  and  the 
experimental  data,  while  step  6 insures  that,  in  the  temperature  range  of  solder  melting,  while  the 
link  passes  through  the  (relatively  small)  temperature  increment  ATmelt  immediately  prior  to  TpuSE, 
the  effective  specific  heat  of  the  link  is  significantly  increased  to  the  value  that  provides  the  amount 
of  latent  heat  to  melt  the  exact  amount  of  solder  required  to  fuse  the  link. 

It  is  of  interest  to  note  that  the  use  of  a relatively-large  effective  value  of  link  specific  heat,  in  step 
6,  for  simulating  the  effect  of  melting  of  the  link’s  solder  on  link  thermal  response  is  analogous  to 
the  use  in  References  [19]  and  [20]  of  a similarly  large  effective  specific  heat  for  simulating,  e.g., 
the  effect  of  evaporation  of  free  water  over  a small  temperature  interval  on  the  thermal  response  of 
fire-exposed  gypsum  wall  board. 

The  above  approach  assumes  that  in  any  relevant  scenario,  once  T,(t)  exceeds  Tmelt  it  continues  to 
increase  monotonically  with  t.  Additional  considerations,  beyond  the  scope  of  this  present  work,  can 
lead  to  modification  to  the  final  model  equations  so  that  they  can  be  used  to  simulate  scenarios  where 
the  link  is:  heated  to  the  point  of  sustained  melting,  with  loss  of  some,  but  not  all  of  its  solder; 
cooled;  reheated  to  the  melt  condition;  etc. 

Using  the  plunge  test  data  for  a final  determination  of  the  system  parameters  for  the 
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configurations  of  Table  2.  The  Seven-Step  Procedure  was  carried  out  for  the  plunge  test  data  and 
corresponding  system  configurations  of  test  1 - 4,  with  TpuSE  = 74  °C  and  ATmelt=  2 CC.  The  results 
are  presented  in  Table  4.  As  indicated  in  the  table,  plots  of  the  final  model  simulation  are  presented 
on  Figures  18-21,  corresponding  to  tests  1 - 4,  respectively. 

As  can  be  seen  in  Figures  1 8 - 2 1 , the  solutions  to  the  final  model  equations,  with  system  parameters 
determined  according  to  the  Seven-Step  Procedure,  yield  simulations  that  provide  excellent 
matches  with  the  data  of  tests  1-  4.  Note  that  in  the  simulations,  when  T1?  the  link  temperature, 
reaches  Tp^,  the  vent  is  predicted  to  open  automatically,  and  no  physical  significance  is  attached 
to  subsequent  (relatively-rapid  rates-of-rise  of)  predicted  values  of  T2. 

Although  the  T,(t)  responses  of  the  essentially  identical  link-mount  configurations  in  tests  1 and  2 
are  closely  reproduced,  and  this  in  spite  of  the  different  flow  orientations,  the  solution  system 
parameters  for  these  test  are  significantly  different,  e.g.,  in  Table  4,  compare  the  RTI/s  (93.71  and 
74.96),  the  RTF’s  (0.2581  and  2.057),  and  the  RTImelt’s  (1681  and  671.4).  Nevertheless,  although 
there  is  a discrepancy  in  the  predicted  critical  value  of  tpuSE,  the  solution  of  the  model  equations  for 
the  set  of  system  parameters  derived  from  either  one  of  the  tests  provides  good  matches  to  the 
experimental  data  of  the  other.  This  is  illustrated  in  Figure  22,  where  the  experimental  data  and  the 
model  simulations  of  both  tests  1 and  2 are  all  plotted  together. 

An  Example  Application  and  Validation  of  the  Model  Equations 


Measured  response  of  the  three-element  link-mount  in  a real  fire  scenario.  This  section  uses 
the  model  equations,  developed  to  characterize  the  three-element  link-mount  system  of  tests  1 and 
2,  to  simulate  the  measured  response  of  this  same  link-mount  system,  used  in  the  vent  design  of  one 
test  fire  scenario  reported  in  [10]. 


Refer  to  Figure  23.  The  fire  scenario  involved  a 1.01  m x 1.01  m square  heptane  spray  burner 
located  0.61  m above  the  floor  in  a facility  with  a flat  square  (30.5  m x 30.5  m)  unconfined  ceiling. 
The  facility  was  identical  to  the  one  used  in  the  tests  of  References  [7]  [8]  [9].  There  it  was  indicated 
that  the  lower  exposed  ceiling  surface  was  made  up  of  0.0159  m-thick  insulating  tiles  with  specific 
heat  753  J/(kg  K),  thermal  conductivity  0.061 1 W/(m  K),  and  density  313  kg/m3.  The  ceiling-to- 
floor  spacing  was  8.2  m.  The  vents  were  1.22  m x 2.44  m x 0.30  m (depth).  They  were  mounted 
in  the  ceiling  with  centers  spaced  at  6. 10  m and  symmetrically  placed  relative  to  the  vertical  axis  of 
the  burner,  i.e.,  the  radius  from  the  burner  axis  to  the  center  of  a vent  was  3.05  m.  Relative  to  the 
axis  of  the  fire  and  the  radial  flow  in  the  ceiling  jet,  the  orientation  of  the  link-mount  system  was 
as  in  Figure  2,  with  expected  parallel  system/flow  orientation  as  in  test  2 (refer  to  Table  2). 


The  flow  of  fuel  to  the  burner  was  controlled  to  yield  a total  fire  energy  release  rate  of 


Q 


FIRE  — 


105[t/(75  s)]  kW,  0 < t < 75  s 
105kW,  75s  <t 


(29) 
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Table  4.  System  parameters  determined  according  to  the  Seven-Step  Procedure  for  the 
system  configurations  and  the  data  of  tests  1-4  (see  Table  2)  with  TFUSE  = 74  °C 
and  ATmelt  = 2 C. 


fit  from  t = 0 to  t = tFUSE 
test  number: 


1 

2 

3 

4 

RTIj/(ms)1/2  = 

93.71 

74.96 

94.56 

71.62 

RTI2/(ms)1/2  = 

160.0 

76.64 

19.26 

NA 

RTI3/(ms)1/2= 

1279 

1507 

2377 

NA 

C12/(ms)1/2  = 

0.2581 

2.057 

2.375 

NA 

C23/(ms)1/2  = 

7.732 

12.30 

4.434 

NA 

Total 

Residuals/N  = 

0.0385 

0.0383 

0.0787 

0.0359 

RTIMEL1/(ms)ia  = 

1681 

671.4 

840.0 

1457 

Figure  no. 

18 

19 

20 

21 

The  fire  was  extinguished  at  150  s. 

In  the  experiment,  the  link  of  the  west  vent  fused  and  the  vent  opened  at  t = 90  s.  During  the  test, 
the  east  vent  link  did  not  fuse  and  the  vent  did  not  open.  (As  indicated  by  thermocouple  data  near 
the  vents,  the  failure  of  the  east  vent  to  open  can  be  explained  by  the  fact  that,  upon  opening  of  the 
west  vent,  the  fire  plume  tilted  immediately  toward  the  west.  As  a result  of  this,  the  ceiling  jet 
velocity  and  temperature  in  the  vicinity  of  the  east  vent  were  reduced  from  their  peak  values  and  the 
temperature  of  its  link  never  achieved  Tp^. 

Simulated  response  of  the  links.  To  use  the  model  equations  to  simulate  the  response  of  the  links 
in  the  above  experiment,  it  was  necessary  to  estimate  u(t),  TG(t),  and  TCEIL(t).  This  was  accomplished 
with  the  use  of  the  fire  model  LAVENT  [6][13][14]. 

The  three-element  link  mount,  including  the  link  itself  is  within  the  vent  cavity,  the  depth  of  which 
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is  0.30  m.  During  the  fire  scenario,  as  the  turbulent  ceiling  jet  flows  past  the  vent  cavity  it  would 
“fill”  the  cavity  with  a turbulent  flow  having  the  same  local  characteristic  velocities  and 
temperatures.  It  is  expected  that  the  characteristics  of  the  cavity  flow  and,  in  particular,  the  values 
of  u(t)  and  TG(t)  local  to  the  link,  are  related  to  the  velocity  and  temperature,  respectively,  of  an 
undisturbed  ceiling  jet  that  would  exist  in  the  absence  of  vent  deployment.  Here  it  is  assumed  that 
u(t)  and  Tc(t)  are  equal  to  the  velocity  and  temperature,  respectively,  of  this  undisturbed  ceiling  jet, 
averaged  from  the  ceiling  surface  and  down  to  the  depth  equal  to  the  depth  of  the  vent,  0.30  m in  the 
present  scenario. 

From  the  above  considerations  and  test  parameters,  the  specified  time  dependent  conditions  of  Eq. 
(10),  necessary  to  solve  the  model  equations,  Eqs.  (8"),  (9),  (12),  and  (27)  were  determined.  These 
are  plotted  in  Figure  24. 

The  model  equations  were  solved  using  the  system  parameters,  as  listed  in  Table  4,  determined  from 
both  tests  1 (normal  flow)  and  2 (parallel  flow).  The  solutions  resulted  in  tpusE  estimates  of  72.5  s 
and  81.0  s.  These  are  to  be  compared  to  the  experimentally  measured  result  of  90  s.  As  mentioned, 
of  the  two  system  characterizations,  the  one  from  test  2,  for  parallel  flow,  which  led  to  the  tp^  = 
72.5  s estimate,  is  expected  to  correspond  more  closely  to  the  system  as  deployed  in  the  particular 
test  under  discussion. 

The  differences  between  measured  and  predicted  values  of  tp^  can  be  attributed  mainly  to 
differences  between  the  actual  vs  expected  history  of  QnRE  as  specified  in  Eq.  (29).  Near-ceiling 
thermocouple  data  reported  in  [10]  for  the  present  test,  and  for  other  tests  with  QHRE  specified  as  in 
Eq.  (29),  indicate  that  in  the  early  several  tens  of  seconds  of  these  tests,  it  is  reasonable  to  expect 
time  shifts  in  QFIRE(t),  i.e.,  delays  of  the  order  of  10-20  s (see  Figures  9 - 11  of  [10]),  leading  to 
corresponding  time  delays  in  the  expected  response  of  link/vent  activation. 

Summary  and  Conclusions 

Model  equations,  suitable  for  general  use  in  compartment  fire  models,  were  developed  [Eqs.  (8")> 
(9),  (12),  (10),  and  (27)]  to  simulate  the  thermal  response  to  arbitrary  fire  environments  of  fusible- 
link  actuated  automatic  fire  vents.  The  method  of  analysis,  which  focused  on  a prototype  three- 
element  link-mount  design,  can  be  extended  to  arbitrary  multiple-element  link  mounts.  Also,  the 
equations  for  the  prototype  three-element  link  mount  were  shown  to  include,  and  be  directly 
applicable  to  the  problem  of  simulating  the  response  of  single-element  simply-supported  links. 

A method  was  developed  to  determine  the  values  of  the  set  of  parameters  that  characterize  a 
particular  multiple  link-mount  design;  five  parameters  in  the  case  of  the  above  three-element  link 
mount.  This  involves:  1)  the  measured  time -dependent  thermal  response  of  the  link-mount  design 
to  exposure  in  a plunge  test,  the  type  of  test  used  to  characterize  the  thermal  response  of  sprinkler 
links;  and  2)  an  analytic  means  of  determining  values  of  the  design  parameters  that  yield  an  optimum 
fit  between  a solution  to  the  model  equations  and  the  temperature  data.  The  method  was  presented 
in  an  earlier  section  of  this  work  entitled  “A  Seven-Step  Procedure  for  determining  system 
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parameters  and  establishing  final  model  equations.” 

In  developing  the  model  equations  within  the  context  of  plunge  test  results,  it  was  determined  from 
the  measured  temperature-time  responses  of  the  links  that,  for  the  particular  link  design  used  in  the 
present  tests,  and,  no  doubt,  for  other  relatively-large-solder-mass  link  designs,  the  time  interval  over 
which  link-solder  melting  takes  place  (i.e.,  the  interval,  prior  to  fusing,  during  which  time  the  link 
temperture  remains  close  to  its  fuse  temperature)  is  a significant  fraction  of  the  total  time  to  fuse. 
For  this  reason,  and  in  contrast  to  the  traditional  way  of  modeling  the  response  of  fusible  links,  it  was 
concluded  that  model  equations  must  account  for  the  solder-melt  phenomenon.  Such  an  accounting 
is  included  in  the  Seven-Step  Procedure. 

The  Seven-Step  Procedure  was  carried  out  successfully  for  four  plunge  tests  involving  three  link- 
mount  systems  and  one  or  two  system/flow  orientations.  In  each  case,  the  determined  values  of  the 
system  parameters  with  the  model  equations  lead  to  excellent  simulations  of  the  plunge  test  data. 
The  resulting  sytem  parameters  and  model  equations  corresponding  to  the  three-element  link  mount 
design  of  the  fusible-link  vent  test  of  Reference  [10]  were  used  successfully  to  simulate  time  of  link 
fusing  and  vent  activation  in  a large-scale  heptane-spray-bumer  fire  environment. 
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Nomenclature 


AcONV.i 

area  of  element  i for 
convective  heat  transfer 

Cij 

Eq.  (14) 

C's 

coefficient  for  conductive  heat 
transfer  from  elements  i to  j 

Cj 

specific  heat  of  element  i 

CMELT 

effective  value  ofjC  during 
solder  melt  phase 

^RAD.i 

Eq.  (12) 

EraD.i.PLUNGE 

Fradi  in  a plunge  test,  Eq.  (18) 

G 

Eq.  (20) 

hi 

average  heat  transfer 

coefficient  of  element  i 

h.AJR 

thermal  conductivity  of  air 

EGONV,i 

characteristic  length  of  element 

m, 

mass  of  element  i 

NF 

number  of  times  data  acquired 
in  plunge  test,  up  to  tpuSE 

Nm 

number  of  times  data  acquired 
in  plunge  test,  up  to  tMELT 

Nu, 

Nusselt  number  of  element  i 

Qfire 

energy  release  rate  of  the  fire 

QraD.i 

Eq.  (9) 

Qradfire,i 

Eq.  (9) 

QcONV.i 

flux  of  convective  heat 
transfer  to  element  i 

QRAD.i 

flux  of  radiative  heat  transfer 
to  element  i 

9radfire 

flux  of  radiative  heat  transfer 
from  the  fire  to  a target 

^RADFIRE.i 

flux  of  radiative  heat  transfer 
from  the  fire  to  element  i 

R 

radius  from  fire  to  target 

Re, 

Reynolds  number  of  of  flow 
around  element  i 

RTI 

response  time  index 

RTI, 

RTI  of  element  i 

RTImelt 

effective  value  of  RTIj  during 
solder  melt  phase 

r 

radial  distance  along  ceiling  to 
link-mount  assembly 

Tame 

ambient  temperature 

^base 

temperature  of  the  base  of  a 
sprinkler 

Tceil 

temperature  of  ceiling 

Texp.! 

experimentally  measured  value 
of  Tj 

Tfuse 

fuse  temperature  of  the  link 

tg 

gas  temperature 

tg. 

Eq.  (3) 

Tlemk 

link  temperature 

Tmelt 

Elink  when  a plateau  in  Tu^is 
initiated 

TwaLL 

wall  temperature  in  plunge  test 

T, 

temperature  of  element  i 

t 

time 

h 

t when  data  is  acquired,  j = 1, 
NF 

^FUSE 

t when  Tj  = 

^MELT 

t w'hen  a plateau  in  Tunk 
is  initiated 

U 

gas  velocity 

UCR 

a critical  value  of  u 

atmelt 

Eq.  (25) 

e 

emissivity 

€i 

e of  element  i 

eCEIL’  €FLOOR 

€ of  element  i,  ceiling,  floor 

eWALL 

e of  wall  in  plunge  test 

A, 

Eq.  (19) 

^RAD 

fraction  of  radiated  to  the 

far  field 

VAIR 

kinematic  viscosity  of  air 

a 

Stefan  Boltzman  constant 

Eq.  (11) 

subscripts 

i 

index  of  element  i;  i = 1,  2,  3 
refers  to  the  link,  tab,  and  bar, 
respectively 
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Figure  1.  The  simply-supported  fusible  link. 
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SPRING-LOADED  (VIA  FUSIBLE  LINK)  VENT 


Figure  2.  A typical  fire  scenario  in  the  NFPRF  program  [9],  including  a sketch  of  the 
vents  with  the  three-element  link  mount. 


THREE-ELEMENT  FUSIBLE  LINK 
(DESIGN  USED  IN  UL  TESTS) 
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The  three-element  link-mount  assembly. 


Figure  4.  Components  of  heat  transfer  to  the  three  elements  of  the  link-mount  assembly 
of  Figure  3. 
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Figure  5.  The  three-element  link-mount  assembly  in  the  plunge  test  of  test  1. 
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Figure  6.  Temperature  data  for  the  link-mount  configuration  and  flow  orientation  of  test 
1 (see  Table  2). 
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Figure  7.  Temperature  data  for  the  link-mount  configuration  and  flow  orientation  of  test 
2 (see  Table  2). 
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Figure  8.  Temperature  data  for  the  link-mount  configuration  and  flow  orientation  of  test 
3 (see  Table  2). 
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Figure  9.  Temperature  data  for  the  link-mount  configuration  and  flow  orientation  of  test 
4 (see  Table  2). 
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Figure  10.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
1 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  t,^^. 
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Figure  11.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
2 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  tMELX. 
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Figure  12.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
3 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  tMELX. 
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Figure  13.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
4 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  tMELX. 
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Figure  14.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
1 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  t^g. 
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Figure  15.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
2 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  t^gg. 
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Figure  16. 


Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
3 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  tFUSE. 
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Figure  17.  Model  simulation  for  the  link-mount  configuration  and  flow  orientation  of  test 
4 (see  Tables  2 and  3),  based  on  optimum  fit  to  data  from  t = 0 to  tFUSE. 
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Figure  18. 


Model  simulation  based  on  the  Seven-Step  Procedure  for  the  link-mount 
configuration  and  flow  orientation  of  test  1 (see  Tables  2 and  4). 


45 


Figure  19.  Model  simulation  based  on  the  Seven-Step  Procedure  for  the  link-mount 
configuration  and  flow  orientation  of  test  2 (see  Tables  2 and  4). 
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Figure  20.  Model  simulation  based  on  the  Seven-Step  Procedure  for  the  link-mount 
configuration  and  flow  orientation  of  test  3 (see  Tables  2 and  4). 
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Figure  21.  Model  simulation  based  on  the  Seven-Step  Procedure  for  the  link-mount 
configuration  and  flow  orientation  of  test  4 (see  Tables  2 and  4). 
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Figure  22.  Model  simulations  based  on  the  Seven-Step  Procedure  for  the  link-mount 
configurations  and  flow  orientations  of  tests  1 and  2 (see  Tables  2 and  4). 
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Figure  23. 


ceiling  (8.2  m above  floor) 

Sketch  of  fire  scenario  for  the  large-scale  two-vent  test  of  [10]. 
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Figure  24.  Simulation  of  TG(t),  u(t),  and  TCEIL(t)  during  the  test  of  [10]  local  to  the  three- 
element  vent  link-mounts. 
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